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Abstract 

We use a one-scale similarity analysis which gives specific relations between the 
velocity, amplitude and width of localized solutions of nonlinear differential equa- 
tions, whose exact solutions are generally difficult to obtain. We also introduce 
kink-antikink compact solutions for the nonlinear-nonlinear dispersion K(2,2) equa- 
tion, and we construct a basis of scaling functions similar with those used in the 
multiresolution analysis. These approaches are useful in describing nonlinear struc- 
tures and patterns, as well as in the derivation of the time evolution of initial data 
for nonlinear equations with finite wavelength soliton solutions. 



PACS numbers: 42.30. Sy, 02.30.Px, 43.35.+d, 47.20.Ky, 05.45.-a, 



1 Introduction 

The nonlinear partial differential equations (NPDE) of physical interest can describe a 
variety of patterns, particle-like traveling solutions, solitons and breather modes in nu- 
clear and particle physics, nonlinear molecular and solid state physics phenomena, and 
features found in nonlinear optics Q. Their solutions are usually localized and demon- 
strate stability in time and in collisions with each other. In the asymptotic domain these 
solutions consist of isolated traveling pulses that are free of interactions and have a shape 
related to the velocity, thus making nonlinear patterns distinct from linear results. In the 
scattering domain, the nonlinear solutions obey some nonlinear superposition principle. 
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Nonlinear dynamics require NPDE which display very strong interaction between the ini- 
tial conditions and the dynamics and involve multiple scales [fj, |3| , being able to produce 
self-similar or fractal-like patterns. Recent examples show that the traditional nonlinear 
tools (inverse scattering, group symmetry, functional transforms) are not always appli- 
cable 0. From the experimental point of view one knows that such patterns generally 
have finite space-time extension and a multi-scale structure. Since the traditional solitons 
or the soliton-like solutions have infinite extent, one needs rather appropriate compact 
supported basis functions to investigate such structures. 

Multi-resolution analysis (MRA), ||, could be a useful method for the construction 
of such nonlinear bases, since the linear harmonic analysis is inadequate for describing 
nonlinear systems. The MRA uses wavelets, which are functions that have a space- 
dependent scale which renders them an invaluable tool for analyzing multi-scale phenom- 
ena. Wavelets have been used in signal processing, problems involving singular potentials 
in quantum mechanics, in discussions concerning q-algebras, and even in nuclear struc- 
ture studies 0. It follows that the use of MRA in the study of NPDE is natural because 
wavelets can analyse nonlinear features like strong variations or singularities. The NPDE 
describe a variety of patterns J2|, features in quantum optics [0, molecular and solid state 
physics phenomena |§ and solitons in nuclear and particle physics || [10 . 



In this paper, wavelet-inspired approaches for localized solutions of NPDE are ex- 
plored. We propose two different formalisms for the scale analysis, and the classification 
of soliton solutions of NPDE. A first method provides relations between the character- 
istics of such solutions (amplitude, width and velocity) without the need of solving the 
corresponding NPDE. The method uses the multi-resolution analysis || instead of the 
traditional tools like the Fourier integrals or linear harmonic analysis which are inade- 
quate for describing such systems. Moreover, the introduction of wavelet analysis in the 
study of NPDE is somehow natural because it can accommodate behaviour diplaying 
strong variations, even singularities, to a smooth behavior. This scale approach has the 
advantage that it does not need the explicit form of the exact solutions. Hence, it is useful 
especially in situations when such solutions are unknown. In section 2 we provide many 
examples, predictions and applications of this one-scale approximation method (OSA) for 
a large class of NPDE with respect to their localized solutions. The NPDE is mapped 
into an algebraic equation relating the amplitude, width and velocity of such signals. The 
results are succinctly presented in tabular form. In section 2 we show an example of the 
construction of a nonlinear basis for NPDE i.e. the construction of a kink-antikink basis 
for the K(2,2) equation, involving nonlinear dispersion. There are many physical reasons 
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favoring wavelets in the construction of such nonlinear bases. The self-similar character of 
the fission process of fluid drops is an example where the same type of singularity occurs 
in any scale |T1], [T2|]. In section 3, we use also this basis (or frame) for the investigation 
of the time evolution of a given initial data profile for the nonlinear-nonlinear dispersion 
K(2,2) equation. The proofs are introduced in two Appendices. 



2 Scale analysis of NPDE 

In this section we introduce a one-scale analysis (OSA) for the NPDE, in terms of their 
localized traveling solutions. The NPDE which describe physical phenomena are usually 
not susceptible to analytic solutions. Moreover, there are examples [Q, |HJ when the 
mathematical tools like the inverse scattering theory or the transformation group method 
are not applicable. When phenomena of interest have many space-time scales, or the 
scale of the process varies in time, the numerical methods may fail, like in the case of 
propagating discontinuities or shock waves. A simple option is the expansion of solutions 
in a basis of appropriate chosen basis. The Fourier series, which has the advantage 
of orthogonality, cannot discriminate the local behavior of phenomena. Moreover, the 
analysis of the Fourier coefficients of a wave u(x, t) is not sufficient for drawing conclusions 
about the scale of the localized structures. 

The so called OSA analysis, described and applied here for localized traveling solutions 
belonging to any type of NPDE, provides algebraic connections between the width L, 
amplitude A, and the velocity V of the solution, without actually solving the equation. 
The procedure consists in the substitution of all the terms in the NPDE, according to the 
rules: 

i. u t — ► — Vu x substitution for traveling solutions, 

ii. u -> ±A, u x -> ±A/L, u xx -> ±A/L 2 . . . , (1) 

and so forth for higher order of derivatives. This substitution in eq.(l ii) is possible only 
for localized (finite extended support) solutions having at least one local maximum (like 
solitons or Gauss functions), if they exist. The advantage of the substitution follows from 
the fact that it enables one to say something about the hypothetical localized traveling 
solutions even if (or especially when) they are not known or discovered. However, this 
averaging method may fail in the case of strictly monotonic solutions, or solutions with 
singularities. This will be made more explicit in the next paragraph and Appendix 1. 
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Since we are interested in traveling solutions, the first substitution reduces the number 
of variables from 2 to 1 so that we are now dealing with an ordinary differential equation 
instead of a PDE. Then, the second substitution transforms the ordinary differential 
equation into an equation in the parameters describing the amplitude, width and velocity. 
Consequently, the NPDE is mapped into an algebraic equation in A, L and V. 

The proof of the method follows from the expansion of the soliton-like solution u(x, t) = 
u(s) with s = x — Vt, in a Gaussian family of wavelets \P(s) = Ne®( s \ where Q(s) is 
a polynomial and N the normalization constant || [L4[]. If we choose Q = —is — y we 
obtain a very particular wavelet with the support mainly confined in the (—1, 1) interval, 
namely ^f(s) = exp[— is — y]/7r 1//4 . We have the discrete wavelet expansion of u 

u(s) = EE^(2 J 3 - k) = £Q, fe %, fc ( S ), (2) 

j k j,k 

in terms of integer translations (k) of which provide the analysis of localization, and 
in terms if dyadic dilations (2- J ) of which provide the description of different scales. 
The idea of the proof is to choose a point where the expansion in eq.(2) can be well 
approximated by one single scale such that u(s) can be approximated with a sum of 
phases in s. In order to reduce the number of scales, the range of the summation should 
be chosen with an eye to the underlying physics. The exact proof of this approach is 
provided in Appendix 1. 

In the following we apply this procedure to three classes of NPDE, namely: disper- 
sive, diffusive and dispersive-diffusive nonlinear equations, and we present the results in 
tabular form. Before describing different examples of NPDE, we point out that for each 
equation we try to present one example of an exact solution (second column Tables 1,2) in 
comparison with the results of the much simplest OSA approach (third column of first two 
tables and second column of Table 3). We stress that we do not choose only that solution 
which fits the similarity analysis. On the contrary, our approach fits a large number of 
solutions, as can be checked directly. 

Our first example of application of the OSA approach is given by the convective- 
dispersive equations, for which the most celebrated example is provided by the KdV 
equation 

u t + uu x + u xxx = 0. (3) 

In general, the stable solutions of the nonlinear-dispersive equations are dependent of 
the initial conditions, through their conservation laws. Consequently, they can gener- 
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ate a large class of patterns, shaped by the balance between nonlinear interaction and 
dispersion, among which the most interesting examples are solitons, breathers and kinks. 

In Table 1 we present examples of pure dispersive NPDE, identified in the first column 
by the form of the equation, and in the second column by a corresponding traveling 
localized solution, if the analytical form is available. Such exact solutions provide special 
relations between L, A and V, which are given in the third column of Table 1. In the 
last column we introduce the results of OSA, namely the relations between these three 
parameters, provided by eqs.(l ii). The usefulness of the approach may be checked, by 
a quick comparison between the second and the third columns. While the results in the 
second column are possible only when one knows the analytical solutions, the results 
presented in the last column, obtained by the OSA, result directly from the NPDE form, 
without actually solving it. 

The case of the KdV equation, eq(3), is described in the first row of the Table 1. The 
OSA method gives a general expression for L = L(A, V), shown in the first line, first row 
of the last column. In the second line of the last column we show that a specification 
of the relation between V and A, together with the general result given in the first line, 
results in a relation between L and A. For L to be related to A only, it results from the 
third column that the velocity V must be proportional to A. In this case we obtain the 
well-known relation (first row, column two) among the parameters in the soliton solution 
V = 2A. Moreover, OSA method allows V to depend on a higher power of A (V ~ A p , 
p > 1). If such a solution could exist, a lower bound for A will occur. Such solitons would 
have only amplitudes higher than this limit, while solitons with a smaller amplitude than 
this limit move with velocity proportional to A. 

Similar results are obtained for the MKdV equation (second row), except that here 
A needs to be proportional to the square root of V in order to have L a function of A 
only. This prediction is again identical with that in the exact solution (second column). 
Moreover, the same relations remain valid even for the solutions of the "compacton" type 



where L = ir/6k, that is L ~ 1/A, like in the Table 1. 

Next example (third row) is provided by a generalised KdV equation, in which the 
dispersion term is quadratic 







u(x, t) 



VS2k cos[k(x - AkH)} 2 



3(1- f cos[k(x - AkH)] 2 )' 



u t + (u 2 ) x + (u 2 ) 



xxx 



0. 



(4) 
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Eq.(4), known as K(2,2) equation because of the two quadratic terms, admits compact 
supported traveling solutions, named compactons ||, [13], [16|, [17|, In general, the 

compactons are obtained in the form of a power of some trigonometric function defined 
only on its half-period, and zero otherwise, in such a way that the solution is enough 
smooth for the NPDE in discussion. In the above example the square of the solution has 
to be continuous up to its third derivative with respect to x. 

Different from solitons, the compacton width is independent of the amplitude and 
this fact provides the special connection with the wavelet bases. The compactons are 
characterize by a unique scale, and it is this feature that makes it possible to introduce a 
nonlinear basis starting from a unique generic function. For eq. (4) the compacton solution 
is given by 



7] C (X - Vt) 



4V 



-cos 



x - Vt 
4 



(5) 



if \x — Vt\ < 2tt and zero otherwise. Here we notice that the velocity is proportional to 
the amplitude and the width of the wave is independent of the amplitude, L — 4. As a 
field of application we mention that the quadratic dispersion term is characteristic for the 
dynamics of a chain with nonlinear coupling. 

The general compacton solution for eq.(4) is actually a "dilated" version of eq.(5). 
That is, a combination of the first rising half of the squared cos in eq.(5), followed by 
a flat domain of arbitrary length A, and finally followed by the second, descending part 
of eq.(5). Actually, this combination is just a kink compacton joined smoothly with an 
antikink one 



Vkak(x - Vt] A) 



0... 



f cos 2 



-Vt 



fcos 2 



-Vt-X 



2tt < x - Vt < 
Vt < A 

\<x-Vt<\ + 2Tt 



(6) 



0. 



In Fig. 1 we present compactons (upper line), kink-antikink pairs (KAK) described by 
eq.(6), both with the same amplitude and velocity (middle line). Although the second 
derivative of this generalized compacton is discontinuous at its edges, the KAK, eq.(6), 
is still a solution of eq.(4) because the third derivative acts on u 2 , which is a function 
of class C3. Finally, we can construct solutions by placing a compacton on the top of a 
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KAK, as in the bottom line of Fig. 1. Such a solution exists only for a short interval of 
time (X/V), since the two structures have different velocities. The analytic expression of 
the solution is given by 

V(x,t) = V kak(x - Vt; A) + L c {x - V't - 2vr) + ^pj x ( X ~ V ^~ 2fK ) ) (7) 

for < t < (A — 4tt)/(V — V) and zero in the rest. Here x( x ) is the support function, 
equal with 1 for \x\ < 1 and in the rest, and V = 3max{r/ c }/4 + 2V. 

For the K(2,2) compacton eq.(5) fulfills some relations between the parameters: A = 
4V/3 and L = 4 . The relations provided by OSA in the last column of the third 



row, predict such relations, and hence also proove the existence of the compacton. That 
is, for a linear relation between the amplitude and the speed, the half-width is constant 
and does not depend on A. Indeed, if we choose V = ±3 A/2 or V = ±5 A/2, we obtain 
L = 4, like for the compacton. The constant value of L = L = const, is a typical feature 
of the K(2,2) compactons. Moreover, it was found numerically that for any compact 
supported initial data, wider than L , the solution decomposes in time into a series of L 
compactons, Fig. 2. For narrower initial data the numeric solution blows up. There is 
no exact or analytic explanation of this effect, so far. The scale relations can give a hint 
in this situation, too, by using the graph of the relation L = L(V, A) provided by this 
qualitative method. In Fig. 3, L from the third column of the third row is plotted versus 
V, for several values of A (larger values of A translate the curves to the right). Above the 
value L = 4 of the half-width of a stable compacton, wider compact pulses produce an 
intersection for each curve (each A) with the axis L = 4, providing a series of compactons 
of different heights, like in the numerical experiments described in 0, [13|, [17| . Below the 
L = 4 line, all the curves L(V) approach L = 0, towards infinite amplitude, explaining 
the instability of the narrower initial data. 

Another good example of the predictive power of the method is exemplified in the case 
of a general convection-nonlinear dispersion equations, denoted by K(n,m) 

Vt + (V n ) x + (V m )^ = 0. (8) 

Compacton solution for any n ^ m are not known in general, except for some particular 
cases. In this case we find a general relation among the parameters, for any n, m, shown in 
the fourth and fifth rows. These general relations L(A, V) approach the known relations 
for the exact solutions, in the particular cases like n = m (fourth row), n = m = 2 (third 
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row), n = m = 3 (first reference in M) and n = 3, m = 2; n = 2, m = 3 (fifth row). These 
results can be used to predict the behavior of solutions for all values of n, m. 

In the following, we present another example of applications of the OSA approach, 
related to a new type of behavior of nonlinear systems. Traditional solitons move with 
constant speed on a rectilinear path (except for the roton [Tl[] which has a circular trajec- 



tory with constant angular velocity). The speed is usually equal to the amplitude scaled 
with a constant. Higher solitons travel faster and there are no solitons at rest (zero speed 
implies zero amplitude). They can travel in both directions with opposite signs for the 
amplitude. The situation is different in the case of compactons, which allow also station- 
ary solutions. When linear and nonlinear disspersion occur simultaneously, like in the so 
called K(2,l,2) equation 

u t + (u 2 ) x + {u) xxx + e{u 2 ) xxx = 0, 
where e is a control parameter, the OSA yields a dependence of the form 



L=^{±A + e)/(V±A), 

which still provides a constant width if V = ±A+2e. In this case, the speed is proportional 
to the amplitude, but can change its sign even at non-zero amplitude. Solutions with larger 
amplitude than a critical one (A crit = =F2e) move to the right, solutions having the critical 
amplitude are at rest, and solutions smaller than the critical amplitude move to the left. 
This behavior was explored in [ 16| , too. However, such a switching of the speed is not 



necessarily a feature of the nonlinear dispersion. A compacton of amplitude A on the top 
of a infinite-length KAK solution of amplitude S 

u(x, t) = Acos 2 -J— ) + 5 ' ( 9 ) 

is still a solution of the K(2,2) equation, with the velocity given by V = f (w + A 

For A = —25 the solution becomes an anti-compacton moving together with the KAK. 
In the case of a slow-scale time-dependent amplitude the oscillations in amplitude can 
transform into oscillations in the velocity. The key to such a conversion of oscillations 
is the coupling between the traditional nonlinear picture (convection-dispersion-diffusion) 
and the typical Schrodinger terms. 

In Table 2 we present another class of NPDE, namely the dissipative ones. These 
equations generalize the linear wave equation (first row) where there is no typical length 
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of the traveling solutions. The wavelet analysis provides the correct expression for the 
dispersion relation (V = c — > k 2 = uo 2 /c 2 ) with no constraint on either the amplitude A 
or on the width L. 

In the second row we introduce the celebrated Burgers equation which represents 
the simplest model for the convective-dissipative interaction. Dissipative systems are 
to a large extent indifferent to how they were initialized, and follow their own intrinsic 
dynamics. We provide in the second column an analytic solution of the Burgers equation. 
For some special of values of the integration constants (2C < V 2 , D = 0) the solution 
becomes a traveling kink 



By applying the OSA approach to the Burger equation (third column) we obtain the same 
relation between amplitude and half-width, like in the case of the exact solution eq.(10), 
providing the velocity is proportional with the amplitude. 

In the following we apply the OSA approach to investigate a nonlinear Burgers equa- 
tion 



called quasi-linear parabolic equation fll8l , and used to describe the flow of fluids in 
porous media or the transport of thermal energy in plasma. The last term describes the 
volumetric absorption (Bremsstrahlung for 7 = 1/2, synchrotron radiation for 7 in the 
range 1.5 — 2, etc). The second term in eq.(ll) describes the convection process and the 
coefficient (3 ranges from to 1 in the case of difussion in plasma, and further to higher 
values, for unsaturated porous medium in the presence of gravity. The existence and 
stability of waves or patterns is strongly dependent on the coefficients a, /x, c, m, k and 7, 
and at this point the OSA can be useful again since there is no general analytic solution 
for eq.(ll). The result of the OSA approach is presented in the third row of Table 2. The 
typical scale of patterns depends on the parameters in the equations and the amplitude 
of the excitations, in a complicated way. 

However, in order to test OSA again, we found a simple class of exact solutions when 
c = 0, presented in the fourt row in Table 2, and expressed as the inverse of a degenerated 
hypergeometric function. In this expression we have A = [ikV a ~ l ((k — l)a a ) _1 , z = 
(a/V)u m ~ 1 and a = (k — l)/(m — 1). The asymptotic behavior of the left hand side of 
the solution given in fourth row, second column, is described by 




(10) 
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If z approaches +00 the solution increases indefinitely like an exponential. For a > 1 
(strong difussion effects), for even k and for even m, the traveling wave u(x — Vt) has 
a negative singularity towards —00 at x + x = T(a + 1)(— l) a < 0. For k odd there is 
also a singularity at x + x$ > 0. These solutions are not likely to provide viable physical 
results. If k is even and m is odd (the singularity is pushed towards imaginary x), or if 
< a < 1, the singularity is eliminated and the solution becomes semi-bounded, like in 
the particular situations investigated in the article fl8 |. In this case, OSA provides again 



the correct relations, since we obtaine the special behavior of the solution if the velocity 
is proportional to the power m — 1 of the amplitude A. Also, we predict the space scale 
of these semi-compact pulses, namely the length L = v ± ar ^ A m-i ■ 

The OSA analysis can be applied in the case of sine-Gordon equation, fifth row of 
Table 2. The solutions with the velocity proportional with L 2 are characterized through 
the OSA approach by a transcendental equation in A, identical with the equation fulfilled 
by the amplitude A of the exact sine-Gordon soliton. 

In the sixth row, we present the cubic nonlinear Schrodinger equation (NLS3) which 
has a soliton solution. This type of equation is applied in nonlinear optics, elementary 
particle physics or in the polaron model in solid state physics 0, §], |9| • Recently, different 
effects of cluster physics could be explained by using the NLS3 equation. In the sixth 
row of Table 2 we present the NLS3 equation together with its one soliton solution of 
amplitude r] , obtained by the inverse scattering method. In the last column we also 
show the relation between the parameters of a localized solution, obtained by OSA. The 
equation for L(A, V) is more general than that one fulfilled by the soliton, and hence is 
related to more general localized solutions. By choosing the velocity proportional to the 
amplitude, we reobtain the =1/770 ~ 1/V typical relations for the soliton given 

in the second column. If a more general solution of the NLS3 equation describes, for 
example, the dynamics of some cluster states, or the dynamics of hard spheres in a hard 
core potential model 

h 2 

then the L parameter gives an estimation for the wavelength of the wavefunctions, or for 
the correlation length in a Bose model 

n n 



~ — -j^^=^^= for small A. 
>2m(E -V) + aA 2 ^2m(E - V) 

For the general case of a NLS equation of order n (seventh row), where a general analytical 
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solution is unknown, the method predicts a special L = L(A, V) dependence, shown in 
the third column and in Fig. 4. Contrary to third order NLS, where the dependence of 
L with A is monotonous for V =~ ±A (n = 3 in Fig. 4), at higher orders than 3, the 
L(A) function has discontinuities in the first derivative. This wiggle of the function (Fig. 
4, n = 4) holds at a critical width, possibly producing bifurcations in the solutions and 
scales. As a consequence, initial data close to this width can split into doublet (or even 
triplet, for higher order NLS) solutions, with different amplitudes. Such phenomena have 
been put into evidence in several numerical experiments for quintic nonlinear equations 



1% |18|, |1J 

The final example of Table 2 is provided by the Gross-Pitaevski (GP) mean field 



equation, which is used to describe the dilute Bose condensate [p0| . The scalar field (or 
order parameter) governed by this equation was shown to behave in a particle manner, 
too, since it can contain topological deflects, namely dark solitons. The space scale L of 
such solutions is important, for both the theory and experiment, since is related to the 
trap dimensions and to the scattering length. In the last row of the Table 2 we give one 
particular solution of a simplified one-dimensional version of the GP equation |yj 

ot \ 2m m J 

where a is the s-wave scattering length and V ex t is the confining potential. In the 
solution provided in the table, the half-width of the exact nonstationary solution is 
L = — p 2 , where v c is the Landau critical velocity, and p = q(t) is the momen- 

tum associated with the motion of this disturbance. It is easy to check that the OSA 
provides a good match with this exact solution, and also L fits the correlation length 



Iq = ^ m/ A-K?i 2 a. We stress that such estimation of the length is also important in nuclear 
physics where one can explain the fragmentation process as a bosonization in a-particles, 
inside the nucleus. Such systems are coherent if the wavelength associated with the clus- 
ter (the resulting L in the GP equation) is comparable with the distance between the 
a-clusters. 

The more complex the NPDE is, the richer the conclusions of the OSA approach. A 
good example of such an analysis is related to the convective-dissipative-dispersive NPDE, 
example provided by the model equation 

ut + a(u m ) x + b(u k ) xx + c(u n ) xxx = (13) 

Here m, k and n are integers and the corresponding terms are responsible of the nonlinear 
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interaction (convective term), dissipation and dispersion The above equation is re- 
lated to weakly nonlinear phenomena, and it occurs in modeling porous medium, magma, 
interfacial phenomena in fluids (and hence applications to drop physics), etc. General 
solutions are difficult to obtain for such a complex equation. We show in the following 
that most of the essential conclusions for the behavior of its solutions can be obtain, in a 
simple way, by the OSA. This approach maps this equation into 

[amA m - x - V)L 2 - fi^A^L + ^A" 1 ' 1 = 0, (14) 

Table 3, first row. The most symmetric case is obtained when either V = (stationary 
patterns) or V ~ A m ~ x . In this situation the condition to have a monotonous dependence 
of L as a function of A is 2k = m + n which yields a scale structure 



/ok 2 ± \//i 2 /c 4 — 4mn 3 (a — V ) \ 

L = A k - m (- ^ [ ^)^A k - m , (15) 

V 2m(a - Vo) / 

where we put V = mVoA m ~ 1 . The condition 2k = m + n is just the condition obtained in 
[ ll| from a scaling approach. This condition assures the universality of the corresponding 



patterns, and it is the unique case in which L depends on a power of A. In any other 
situation, this dependence is more complicated and introduces singularities which broke 
the self-similarity. In the above cited paper, the author finds out the condition for mass 
invariance at scaling transformations asm = n + 2 = k+l. In our case we just have to 
request the product AL (which gives a measure of the mass, or volume of the pattern, 
like in the case of one- dimensional solitons) to be a constant. This gives the condition 
k—m+1 = which, together with the general invariance condition 2k = m+n, reproduces 
m = n + 2 = k + l. In this case we have patterns characterized by a width 



uk 2 ± Ju 2 k 4 - Amn 3 (a-V ) 

L - \ mA( a-vJ ^ 

If a ~ Vo the width approaches n 3 /A/j,k 2 . In order to make L independen of A, like in 
the compacton case, we need m = k, which together with the first invariance condition 
2k = m + n, yields m = n = k. This is the exceptional case when the dissipative and 
dispersive processes have the same scaling, resulting form the invariance of the eq.(13) 
under the group of scales. Finaly, if we choose L ~ V we obtain the condition k + 1 = 2m 
which (together with 2k = m + n) is the condition for spiral symmetry and occurence of 



similarity structures flql . 
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The above comments are not intended to be a complete study of such a complex 
equation, but the just a proof of how many conclusions one can obtain from the simple 
equation in A, L and V, eq.(14). 

The next example is provided by one of the most generalized KdV equation, which is 
generated from the Lagrangian []13| 



C(n, I, m,p) 



(16) 



where a, (3 and 7 are parameters adjusting the relative strentgh of the interactions, and 
n,l,m,p are integers. For example, for 7 = one re-obtains the K(2,2) equation, and 
for 7 = m = 0,p=l one obtains the KdV equation. The associated Euler-Lagrange 
equation in the function <p x = u(x,t) — ► u(x — Vt) = u(y), reads after one integration 

Vu = -^—u p+1 - (3mu m ~\u y ) 2 + 2(3{u m u y ) y + ^u n -\u y )\u yy f 

- j(u n (u y y-\u yy ) 2 ) y + l{u n {u y ) l u yy ) yy + C, (17) 

where C is the integration constant. By using the OSA we obtain the following important 
result, expressed in the second row of Table 3: The unique case when such an equation 
allows compact supported traveling solutions is when m = p = n + r,C = and V = 
VqA 171 . This result is in full agrement with the variational calculation in 

Both eqs.(13) and (17) are rather more qualitative than capable of modeling measur- 
able phenomena. That is why we introduce now a more general model equation, in the 
form 

u t + f(u) x + g{u) xx + h(u) xxx = 0, (18) 

where f,g and h are different iable functions of the the function u(x,t) itself. The OSA 
approach gives the equation 

-V + f(A) + A 9"( A )+9'{A) + A 2 h"'(A) + 3Ah"(A) + h'(A) = Q (ig) 

A general analysis of eq.(19) is difficult, and the best ways are numerical investigations 
obtained for particular choices of the three functions. We confine ourselves here only 
to show that the class of solutions which have similarity properties are those for which 
V = V f'(A). In this case eq.(19) can be reduced to 

L 2 f'(l - V Q ) + L(Ag" + g') + A 2 ti" + 3Ah" + h' = 0, (20) 
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case which is presented in the third row of Table 3. This last relation can be used for 
different purposes. For example, given a certain type of dispersion and difusion (g, h 
fixed), we can estimate for what types of nonlinearity (/) the width L will have a given 
dependence with A. Or, if we know for instance f(u) = foU qi and h(u) = h u q2 , we 
can ask what type of diffusion g we need, to have constant scale (width) of the patterns 
(waves), no matter of the magnitude of the amplitude A. In other words, which is the 
compatible diffusion term, for given nonlinearity-dispersion terms, which provides fixed 
scale solutions. The result is obtained by integration eq.(20) with respect to g{u) 

9(u) = -^(l + q 2 + -^-)u- - L/ ° (1 ~ yo) ^ + C 3 Logu + C A , (21) 
L \ q 2 — 1 / qi - 1 

where 6*3,4 are constants of integration. In a similar way one can check the existence of 
different other configurations by solving eq.(20), or more general, eq.(19). 

A last application of this method, occurs if the KdV equation has an additional term 
depending on the square of the curvature 

u t + uu x + u xxx + t(u 2 xx ) x = 0. (22) 

This is the case for extremely sharp surfaces (surface waves in solids or granular materials) 
when the hydrodynamic surface pressure cannot be linearized in curvature. Such a new 
term yields a new type of localized solution fulfilling the relations 



AeA 



\ ±^Jl-8eA(A±V) - 1 

If we look for a constant half- width solution (compacton of 1/L = a) we need a dependence 
of velocity of the form V = (l + a 2 e/8)A + l/8eA + a/4:. There are many new effects in this 
situation. The non-monoton dependence of the speed on A introduces again bifurcations 
of a unique pulse in dublets and triplets. Also, there is an upper bound for the amplitude 
at some critical values of the width. Pulses narrower than this critical width drop to zero. 
Such bumps can exist in pairs of identical amplitude at different widths. They may be 



related with the recent observed "oscillations" in granular materials |L6|, p22| . 

The examples presented in Tables 1-3 prove that the above method provides a reliable 
criterion for finding compact suported solutions. The reason this simple prescription works 
in so many cases follows from the advantages of wavelet analysis on localized solutions. 
We stress that this method has little to do with the traditional similarity (dimensional) 
analysis j|, |T3|, [H| [T7|, In the latter case one obtains relations among powers of 
A, L and V, not relations with numeric coefficients like those found in our method. 
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3 The frame of KAK pairs 



In the following we investigate the posibility of construction of a nonlinear frame (an over 
determined or incomplete basis) by using some compact solutions of the K(2,2) equation. 
The high stability against scattering of the K(2,2) compactons, or compacton generation 
from compact initial data, suggest they may play the role of a nonlinear local basis. We 
know from many numerical experiments |12|, [13|, [TJ|], that any positive compact initial 
data decomposes into a finite series of compactons and anticompactons. This suggests 
that an intrinsic ingredient for a nonlinear basis could be the multiresolution structure of 
the solutions, similar with the structure of scaling functions in wavelet theory. 

The compactons given in eqs.(5,6) have constant half-width and hence describe a 
unique scale, which can cover all the space by integer translations. From the point of 
view of multi-resolution analysis, the K(2,2) equations act like a L-band filter, allowing 
only a particular scale to emerge for any given set of initial condition. To each scale, from 
zero to infinity, we can associate a K(2,2) equation with different coefficients. However, 
the compacton solution is not the unique one with this property. For a given K(2,2) 
equation, we can thus extend the scale from L to any larger scale. These more general 
compact supported solutions are still C 2 (R) and are combinations of piece- wise constant 
and piece-wise cos 2 functions. The simplest shape is given by a half-compacton prolonged 
with a constant level, that is a kink solution. The basis solution is a kink-antikink (KAK) 
compact supported combination, Fig. 1. Such kink-antikink pairs of different length, can 
be associated with other compactons, or KAK pairs, one on the top of the other 



Vcomp+KAK(% ~ Vt\ A) — < 



0... 

w 

3 

W 
3 ' 
4V 

3 :! 
4V 

3 ' 
4 V 



COS 



< x 



x-Vt 



-2ir <x-Vt<0 
Vt<5 



+ UV-2V) cos 2 



x-V't 



5<x-Vt<5 + A-n (23) 



5 + 4-K<x-Vt<\ 



3 
0... 



COS 



x-Vt-X 



\<X-Vt<X + 27T 



where 5 < A characterizes the initial position (at t = 0) of the top compacton, with respect 
to the flat part of the KAK solution. The amplitude 4(V — 2V)/3 of the compacton, and 
the amplitude 4V/3 of the KAK, are related to their velocities V and V, respectively. 
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The length of the flat part, A, is arbitrary. A compound solution is not stable in time since 
the different elements travel with different velocities. The total height of the compacton 
is 4(V — V)/3. Since the higher the amplitude is, the faster the structure travels, the 
top compacton moves faster than the KAK, and at a certain moment it passes the KAK. 
Because the area of the solution is conserving, such a compound structure decomposes 
into compactons and KAK pairs. Similar and even more complicated constructions can 
be imagined, with indefinite number of compactons and KAK's, if one just fulfills the C3 
continuity condition for the square of the total structure. Such structures, defined at the 
initial moment can interpolate any function, playing a similar role with wavelets or spline 
bases. It has been also proved that the KAK solutions are stable, by using both a linear 
stability analysis and Lyapunov stability criteria. 

For a given K(2,2) equation, the compacton solution, eq.(5) and in addition the family 
of KAK solutions, eq.(6) can be organized as a scaling functions system. They act like a 
low-pass filter in terms of space-time scales and give the opportunity to construct frames 
of functions from the wavelet model f|5l H, IHf . 



With the notation from Appendix 2, and from eq.(32), we have the elements of the 
frame 

r)k,j( x ) = Vkak(Tr(x - 2 3 Vt - fc), 2 J - l)| i=0 , 

where t — means that we neglect the time evolution, but the amplitude is still amplified 
with a factor of 2- 7 , in virtue of relation r] max = 4V/3. We can now expand any initial 
data for the K(2,2) equation in this frame 

M x ) = c k,jVkA x )> ( 24 ) 
k j 

and taking into account that 

^1,? [ = o otherwis6; (25) 
we can show that the square of the initial data can be linearized by 

u 2 (x) = Yl 12 c kjC k >,j> 

fej j'>j k'ei 



x 



1 1 



EE- E Va(i 1 ,i 2 ,...,i 1 ,_Aj')Vk', j ', (26) 



il=0i2=0 iji_j=0 



16 



where / is the range of k' described in the first line of eq. (25), and 

i'—i 

i=i 

+ k2 U'-^ U '- j ^- j+i " . (27) 

From this relation and from eq.(25) we notice that in eq.(27) the unique nonzero terms 
are those for which a(ii, ii, ij'-j) = k' with k' e /. Hence the initial data is expanded 
in different scales at different translations. The translations are mutually orthogonal so 
they do not give a contribution to the square. In the multiplication of two different scales 
in the expression of the square, we reduce the wider scale in terms of linear combination 
of the narrower ones, by using the two-scale equation, eq.(31) in Appendix 2. All the 
nonzero terms in this product are of the order (2~- J — l)/(2~ jf — 1) ~ 2 J ~K This number, 
given by the number of solutions of the equations 17(21,22, V-j) — k', with k! G /, is 
much smaller than the initial number of terms, where from the advantage of the frame. 
This is the advantage of treating nonlinear problems with a basis that has a scale criterion. 

Another application of the KAK basis occurs when one needs to understand the dy- 
namics of the initial data for the K(2,2) equation. Many numerical simulations show that 
in the case when the width of the initial data is larger that L cornpacton , the initial shape 
decomposes into a finite number of compactons having the same width and different am- 



plitudes |4], [L3|, [L7| p!8| . If we take such an initial pulse, Fig. 5a, and we want to see its 
time evolution untill the break up process, we have first to expand it into the KAK frame, 
Fig. 5a. Then we just let the system develope and the different KAK's move with their 
correspon velocity, eq.(23), providing the new shapes, Fig. 5b and 5c, at different mo- 
ments. We notice a good agreement between these theoretical calculation and the above 
quoted numerical experiments. At this stage we can not yet predict the way that the 
KAK breaks back into compactons. However, we notice that the area of such structures 
in the phase space is the Poincare invariant and should give a hint towards the breaking 
process. This invariant is nothing that the u 3 (denoted D2 invariant in [4]) invariant of 
the K(2,2) equation. 

4 Comments and conclusions 

First, we make a general statement concerning the compact solutions of one-dimensional 
NPDE. A one-dimensional dynamical model is described by a general NPDE equation 
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dtu = 0(x, d x )u where O is a nonlinear differential operator. By taking into account only 
traveling solutions, this NPDE reduces to a nonlinear ordinary differential equation in 
the coordinate y = x — Vt for an arbitrary velocity V. If u(y) is a compact supported 
solution it follows that it is not unique with respect to fixed initial compact data. Indeed, 
if we choose the initial data such that the function and its derivatives (up to the requested 
order) are zero in the neighborhood of a certain point y$ of the y axis, these conditions can 
be fulfilled by any translated version of a compact supported particular solution, placed 
everywhere on the axis outside this neighborhood. Consequently, for such initial data, 
the solution is not unique. This result shows that the compact supported property of the 
initial data and of the solution implies its non-uniqueness. 

Since we can transform the NODE into a nonlinear differential system of order one 

U y = F(y,U), U=(u,u y ,...), (28) 

we can apply the fundamental theorem of existence and uniqueness to solutions of eq.(28), 
for given initial data U(yo) = Uq. If the function F in eq.(28) fulfills the Lipschitz condition 



(its relative variation is bounded) than, for any initial condition, the solution is unique [2| 
Since any linear function is analytic and hence Lipschitz, we conclude that only nonlinear 
functions F allow the existence of compact supported solutions. Thus, a compact soliton 
implies non-uniqueness in the underlying NPDE, which implies non-Lipschitzian structure 
of the NPDE and hence the existence of nonlinear terms. 

In this paper we introduce new physical applications for wavelets, that is the study of 
localized solutions of nonlinear partial differential equations. The existence of compactons 
underlines a common feature of NPDE, discrete wavelets, and also finite differences equa- 
tions. We propose a new scale approach for the similarity analysis and clasiffication of 
soliton solutions, without the need of solving the corresponding NPDE. Also, we proved 
that starting from any unique soliton solution of a NPDE, we can construct a frame of solu- 
tions organized under a multiresolution criterium. This approach provides the possibility 
of constructing a nonlinear basis for NPDE. We show that frames of self-similar functions 
are related with solitons with compact support. In addition, we notice the evidence that 
compactons fulfil both characteristics of solitons and wavelets, suggesting possible new 
applications. Such a unifying direction between nonlinearity and self-similarity, can bring 
new applications of wavelets in cluster formation, at any scale, from supernovae through 
fluid dynamics to atomic and nuclear systems. The scale approach can be applied with 
success to the physics of droplets, bubbles, patterns, fragmentation, fission and fusion. 
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5 Appendices 
5.1 Appendix 1 



We base our proof on the similar proposition in p5| , excepting that here we use Gaussian 
filtering instead of Morlet. We start form the the discrete wavelet expansion of the signal 
u(s) given in eq.(2) 

j k j,k 

in terms of integer translations (k) and dyadic dilations (2- 7 ) of the \I/ wavelet. In order 
to reduce the number of scales needed, the range of the summation should be choosen as 
a function of s. We use in the following the asymptotic formula describing the pointwise 
behavior of the Gaussian wavelet series around a point Sq of interest [2fJ. For a chosen 



So and scale j, there is only one k and |e| < 1 such that the support of the corresponding 
^j t k contains this point, k = 2 J so + e. We can express the solution and its derivatives in 
a neighborhood of this point 

u(s ) « *HE cf j,J*«+f = H M j( s o), 

3 3 

u x (s ) « -i*(-e)^2^C i>2Jso+e = -i'E&u&o), (29) 

3 j 

for e being chosen enough large compared to one, and where Uj(so) = ^(—e)Cj i2 3so+e ~ 
^(0)Cj 2 j so ... Since the coefficient 1/2- 5 represents the scale for each Gaussian wavelet, 
we can define it as a characteristic half- width Lj. Also, we finally have for the n-th order 
derivative in s 

u xx ... x (s )*j:^-. (30) 

3 L 3 

Eq.(30) is the multi-scale generalization of the simpler formula ii in eq(l). With eq.(30) in 
hand we can investigate the structure of hypothetic soliton solutions of NPDE, by choosing 
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So in the neighborhood of the maximum value of the solution, u(sq) = A. Around this 
maximum, such solutions can be described very well by a unique scale L, and hence the 
solution and its derivatives can be approximated with the corresponding dominant term, 
by the substitutions in eq.(l). 



5.2 Appendix 2 

For the sake of simplicity we will renormalize the coefficients of the K(2,2) equation 
such that the support of the simple compacton is one. That is, we take r] c (x,t) = 
r}k a k(7i(x — Vt),0) on the interval \x — Vt\ in [—1/2,1/2]. We construct a multireso- 
lution approximation of L 2 (R), that is an increasing sequence of closed subspaces Vj, 
j G Z, of L 2 (R) with the following properties H, Tl 



1. The Vj subspaces are all disjoint and their union is dense in L 2 (R). 

2. For any function / G L 2 (FL) and for any integer j we have f(x) G Vj if and only if 
D^ 1 f(x) G Vj-i where D^ 1 is an operator that will be defined later. 

3. For any function / G L 2 (Tl) and for any integer k, we have f(x) G Vq is equivalent 
to f(x - k) G Vq. 

4. There is a function g(x) G Vo such that the sequence g(x — k) with k G Z is a Riesz 
basis of Vq. 

In the case of compact solutions of K(2,2) of unit length, we chose for the space Vq that 
which is generated by all translation of r] c with any integer k. The subspaces Vj for 
j > are generated by all integer translations of the compressed version of this function, 
namely, by r\ kak (^Pi\(x — Vt), 0). The subspaces Vj for j < are generated by all integer 
translations of the KAK solution of length A2 J — 1. For example, V-\ is generated by 
Vkaki^ix — 2 J Vt),0). The spaces Vj,j > are all solutions of K(2,2); the others are not. 
The function g(x) is taken to be fjkaki^ix — Vt), 0). It is not difficult to prove that these 
definitions fulfill restrictions one, three, and four. As for the second criterion, we define 
the action of the operator D^ 1 f{x) = f(2x) if f(x) G Vj with a j positive integer, and 
D^a^ix - 2'Vt),2-i - 1) = Vkak(^(x - 2-1 +1 Vt),2-* +1 - 1) for negative j. In 
conclusion, we construct a frame of functions made of contractions of compactons and 
sequences of KAK solutions. We can write the corresponding two-scale equation which 
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connects the subspaces 

Vkak(^(x - Vt), 1) = i] kak (7r(x - Vt), 0) + i] kak (7r(x -Vt-1),0). (31) 

We will denote generically by rjkj the elements of this frame, that is 

Vkj(x) = VkakMx ~ Wt - fc), 2 j - 1)U =0 , (32) 

where t — means that we neglect the time evolution, but the amplitude is still amplified 
with a factor of 2 J , in virtue of relation i] max = 4V/3. 



21 



References 



[1] E. Alfinito, et al, eds., Nonlinear Physics, Theory and Experiment (World Scientific, 
Singapore, 1996); M. Remoissenet, Waves called Solitons ( Springer- Verlag, Berlin, 
1999). 

[2] L. P. Kadanoff, From Order to Chaos (World Scientific, Singapore, 1993). 

[3] J. U. Brackbill and B. I. Cohen, eds., Multiple Time Scales (Academic Press, Inc., 
Orlando, 1985); C. Etrich, U. Peschel, F. Lederer and B. A. Malomed, Phys. Rev. E 
55 (1997) 6155; 

[4] P. Rosenau and J. M. Hyman, Phys. Rev. Let. 70 (1993) 564; B. Dey, Phys. Rev. E 
57 (1998) 4733. 

[5] I. Doubechies and A. Grossmann, J. Math. Phys. 21 (1980) 2080; A. Grossman and 
J. Morlet, SI AM J. Math. Anal. 15 (1984) 72; 

[6] G. Kaiser, A Friendly Guide to Wavelets (Boston, Birkhauser, 1994); A. Ludu and J. 
P. Draayer, in Group 22. Proceedings of The XXII Int. Colloq. on Group Th. Methods 
in Phys., Eds. S. P. Corney, R. Delbourgo and P. D. Jarvis (Int. Press, Cambridge, 
Massachusetts, 1999). 

[7] J. R. Taylor, ed. Optical Solitons - Theory and Experiment (Cambridge University 
Press, Cambridge, 1992). 

[8] A. S. Davydov, Solitons in Molecular Systems (Reidel, Dordrecht, 1990). 

[9] A. Chodos, ed., Solitons in Nuclear and Elementary Particle Physics (World Scien- 
tific, Singapore, 1984); A. Ludu, A. Sandulescu and W. Greiner, J. Phys. G: Nucl. 
Part. 21 (1995) L41; V. Kartavenko and W. Greiner, Int. J. Mod. Phys. E, 7 (1998) 
287. 

[10] R. Rajaraman, Solitons and Instantons (North-Holland, Amsterdam, 1984) 

[11] A. Ludu and J. P. Draayer, Phys. Rev. Lett. 10 (1998) 2125; J. M. Lina and M. 
Mayrand, Phys. Rev. E 48 (1994) R4160. 



22 



[12] J. L. Bona, et al, Contemp. Math. 200 (1996) 17 and Phil. Trans. R. Soc. Lond. A 351 

(1995) 107; G. Zimmermann, Proceedings Int. Conf. on Group Theoretical Methods 
in Physics, G22 (Hobart, 10-18 July, 1998) in press and private communication. 

[13] F. Cooper, J. M. Hyman and A. Khare, Compacton Solutions in a Class of Gener- 
alized Fifth Order KdV Equations in press. 

[14] C. K. Chui and A. Cohen, in Approximation Theory VII, E. W. Cheuey, C. K. Chui 
and L. L. Schumaker, eds. (Academic Press, Boston, 1993); D. Han, Y. S. Kim and 
M. E. Noz, Phys. Lett. A 206 (1996) 299. 

[15] C. N. Kumar and P. K. Panigrahi, Preprint \solv-int/ 9904 QEQ - 

[16] P. Rosenau, Phys. Lett. A 211 (1996) 265 and Phys. Rev. Lett 73 (1994) 1737. 

[17] J. M. Hyman and P. Rosenau, Physica D 123 (1999) 502. 

[18] P. Rosenau, Physica D 123 (1999) 525; 8D (1983) 273. 

[19] V. G. Kartavenko, A. Ludu, A. Sandulescu and W. Greiner, Int. J. Mod. Phys. E 5 

(1996) 329. 

[20] M. H. Anderson et al Science 269 (1995) 198. 

[21] Th. Busch and J. R. Anglin, Phys. Rev. Lett., 84 (2000) 2298. 

[22] P. B. Umbanhowar and H. L. Swinney, Nature 382 (1996) 793; 

[23] V. I. Karpman, Phys. Lett. A 210 (1996) 77; M. I. Weinstein, Comm. Math. Phys. 
87 (1983) 567; A. Khare and F. Cooper, Phys. Rev. E 48 (1993) 4853; B. Dey, C. 
N. Kumar and A. Khare, Phys. Lett. A 223 (1996) 449. 

[24] W. Hurewicz, Lectures on Ordinary Differential Equations (M.I.T. Press, Cambridge, 
Massachusetts, 1985) pp. 5-10. 

[25] A. Ludu, G. Stoitcheva and J. P. Draayer, Int. J. Mod. Phys. E in print. 

[26] N. N. Reyes, J. Approx. Th. 89 (1997) 89. 



23 



Figure Captions 



Fig. 1 

Examples of compactons, kink-antikink pairs, and mixed solutions of the K(2,2) 
equation, together with their velocities. 

Fig. 2 

A finite series of K(2,2) compactons emerging from initial compact data of width 
larger than of the compacton. 

Fig. 3 

The half-width L versus velocity V for the K(2,2) equation, for different amplitudes 
A. Widths larger than L cornpacton = 4 behave different than narrower widths, with 
L < 4. Amplitude increases from left to right, in the range 0.01-0.85. 

Fig. 4 

The half-width L versus amplitude A, for the third (n — 3), forth (n = 4) and fifth 
(n = 5) order NLS equation, in two V = ±A cases. We notice that the higher order 
(n > 3) NLS equations have bifurcations. 

Fig. 5 

The dynamical evolution of the expansion in the KAK basis of a finite initial pulse 
of width L much larger than a compacton, drawn for three moments if time. 
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Table 1: Traveling localized solutions for nonlinear dispersive equations. 



NPDE 


Analytic solution and the 




OSA 




relations among parameters 




approach 


u t + 6uu x + u xxx = 


A sech 2x ~^ t ; L = ^2/A, 


L = 


= \V±QA\- 1 / 2 




V = 2A 


If V 


- A, L ~ A~ 1/2 


u t + u 2 u x + u xxx = 


Asech x -p-, L — 1/A, 


L = 


\V±6A 2 \~ 1 / 2 




A = VV 


If V 


~A 2 ,L~A- 1 



u t +(u 2 ) x + (u 2 ) xxx = Acos 2 ^P, if \{x-Vt)/A\ <tt/2; L = (j^) 

L=4 



2/ x-Vt 



l 

71-1 



if b -Vt\< 

1 1 — n— 1 



n(n 2 +l) 
ct+n. 



1/2 



and else; 

£ — 4n /i 2Vn 



_ 4 ZVn 

1) ' 71 ~~ n+l 



if V = ctA n - x 



u t + {u n ) x + {u"% 
n 7^ m 



unknown 
in general 



n(n 2 +l)A" 
V±mA m - 1 



l\V2 
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Table 2: Traveling localized solutions for nonlinear diffusive equations. 



NPDE 


Analytic solution and the 
relations among parameters 


OSA 
approach 


U X x — ^ u tt = 


/c 2 = cu 2 /c 2 


V = c 
A, L arbitrary 


Ut + uu x — u xx = 


ytv — v tan^vu ~~ ^ — 2 r-^J 


if y ~ A, l ~ 1/4 


ut + a(u m ) x - n{u k ) xx + cu 1 
=0 


only particular cases 
known 


cA^L 2 + (V ± ami" 1 ' 1 )! 
±lik 2 A k - 1 = 


u t + a(u m ) x - u{u K ) xx = 


— Az a ±Fi(a, a + l,z) = x + x 


L = ^ A k m , 

am— a 

if V = aA m - x 


u xt — sin u = 


A tan 7 e £ 


±^4 = sin A 
If V = L 2 ,A = sinA 


itf t + tf xa: + 2|tf| 2 tf = 


■r]Qe i(u}t+kx ' ) sech[i] (x - Vt)]; 
L = I/770 


T ±V±y/\V*-4A*\ 

IiA~V, L = l/A 


^ + ^ + 1^-^ = 


unknown in general 




i% = -|A* 
+ [a|^| 2 + V(:r) - 


ip + \/v 2 - p 2 x 
tanh[a^ 2 -p 2 (x-^))] 


L = {aA 2 ±V- iy 1/2 
If V ~ ±1,L~ I / (Ay/a) 
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Table 3: Traveling localized solutions for dissipative-dispersive equations. 



The NPDE equation 


OSA approach 


Ut + a(u m ) x + b(u k ) xx + c(u n ) xxx = 0; 


j Am-k fJ.k 2 ±y / / J , 2 k 4 ~4mn'- i (a-Vo) 
^ — A 2m(a-V Q ) 

if V = mV A k - 1 


Vu = ^uP +1 - fimu" 1 - 1 ^) 2 + 2f3{u m u y ) y 
+^-\u y )\u yy f - ^{u n {u y ) l -\u yy ) 2 ) y 

+ l{u n {Uy) l Uyy)yy + C 


2L l+4 ((n + l + l)V -a) 
-2L l+2 (l + n + l)(l + n + 2){3 
-(l + n+ 1)(2 + 2n 2 + 31 + I 2 + n(5 + 3Z))7 
if C = 0, V = V A m and m = p = n + I 


u t + f(u) x + g(u) xx + h{u) xxx = 


L = -[<?' + Ag" T [{Ag" + g') 2 

\ 1/2-1 

-4/'(l - V )(A 2 ti" + 3Ah" + h')j 

x(2A 2 h" f + 6Ah" + 2h')- 1 J 
if V = Vof'(A) 
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